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We show that the dynamical behavior of a coupled map 
lattice where the individual maps are Bernoulli shift maps can 
be solved analytically for integer couplings. We calculate the 
invariant density of the system and show that it displays a 
nontrivial spatial behavior. We also introduce and calculate 
a generalized spatiotemporal correlation function. 
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I. INTRODUCTION 

The study of temporal chaos in low-dimensional sys¬ 
tems, some of which can be described by low-dimensional 
maps was extremely beneficial for the understand¬ 
ing of turbulence. In 1984 coupled map lattices were 
introduced into the physical literature as a tool for study¬ 
ing spatiotemporal chaos in spatially extended, i.e. high¬ 
dimensional systems 0. They consist of spatially cou¬ 
pled low-dimensional maps and represent dynamical sys¬ 
tems that are discrete in space and time, but continuous 
in the state variable. They serve as models for coupled 
Josephson junctions, excitable media, population dynam¬ 
ics, neural dynamics and turbulence ||]. Although Buni¬ 
movich and Sinai mathematically proved a number of 
statements regarding the appearance of coherent struc¬ 
tures from spatiotemporal chaos ||, most results in the 
field have been obtained by numerical simulations 
While the study of temporal chaos has greatly profited 
from the existence of simple maps like the Bernoulli shift 
map and the cat map |^],^|, which can be solved explicitly 
(for integer expansion rates), thereby making the mecha¬ 
nisms of mixing and temporal chaos understandable, no 
investigation of this type has been provided up to now 
for the problem of spatiotemporal chaos. Here we present 
a solvable model for spatiotemporal chaos which can be 
solved in the sense of the Bernoulli map. 

We start from the model equations || 

4 +1 = (1 - e)f(zl) + I [f(z* + 1 ) + /(4-i)] , (1) 

where the index i runs over the N sites of a discrete 
lattice, e measures the strength of the spatial coupling 
between neighboring sites and f(zj) is a local map which 
determines the evolution of the continuous variable z\ at 
discrete time-steps t = 0,1,2,-- The extension of eq. 
(|I) to more than nearest neighbor coupling and arbitrary 
dimensions is straightforward. We first write it in the new 
variables x\ = f(z*) as 


.,*+1 


= / 




( 2 ) 


In order to obtain a solvable model we choose f(x) = 
(ax) mod 1, i.e. the Bernoulli shift map with integer 
stretching factor a. To make our approach for the spa¬ 
tially extended system more transparent we will calcu¬ 
late in Section |I| the invariant density and a temporal 
correlation function for the single Bernoulli shift map. 
In section 0 we will show that the coupled map lattice 
is solvable for special values of the coefficients a and e in 
eq. (||) and we will also calculate the invariant density 
and a spatial and a spatiotemporal correlation functions. 
Finally in section IV we discuss our results and indicate 


directions of further research. 


II. PROPERTIES OF A SINGLE BERNOULLI 
SHIFT MAP 

First we recall that the single map x t+1 = (ax*) mod 1 
can be solved as x* = ( a t x°) mod 1 because 

x 2 = {a [(ax°) mod l] } mod 1 

= {a \ax° — k° ] } mod 1 = (a 2 x°) mod 1 (3) 

where k° is an integer which represents the action of the 
modulo function. For the last equality sign in to hold 
we needed the fact that a is an integer such that ak° 
becomes again an integer, which can be dropped within 
the last modulo function. 

Since the modulo function confines the variable x* to 
a circle we could view the Bernoulli shift map as a linear 
map x t+1 = ax* where the variables live on a unit circle, 
i.e. on a 1-torus. We shall see below that we can view our 
coupled map system as a linear map acting on variables 
confined to an IV-torus, where N is the number of lattice 
sites. 

The invariant density p(x) of the simple Bernoulli shift 
map measures the distribution of x values on the attrac¬ 
tor generated by the map and is well known to be a con¬ 
stant ||. We can obtain this result by noting that p(x) 
is defined on an unit circle, i.e. it is periodic in x and 
therefore can be represented as a Fourier series 

p(x)=J2me 2 ” ikx , ( 4 ) 

k 

where k takes only integer values k = 0, ±1, ±2, • • •. The 
invariant density p(x) evolves from an initial distribution 
p°(x) according to the Frobenius-Perron equation [|J 


1 





p* (x) = / dx'5 [x — (a*x*) mod l] p° ( x ') (5) 

Jo 

and is defined as the long-time limit p{x) = lim^oo p*(x). 
In order to solve eq. (|j) we use ([|) and the fact that the 
Bernoulli shift map becomes a linear map on a torus, i.e. 
exp(27ri [( a*x ) mod 1]) = exp(27rm t x), and obtain 

p*(k)=p°(a*k). ( 6 ) 

If we make the reasonable assumption that the initial 
distribution p°(x) is non-singular, then lim^j-oo p°(k) = 
0. This means that in eq. ([|) all Fourier coefficients of 
p*(k) tend to zero in the infinite-time limit, except the 
one, which belongs to k = 0. Since 

p°(0)= fdxp°(x) = l (7) 

Jo 

this yields p{k) = 5k, o and p(x) = 1. 

In a similar fashion we can now define and calculate 
the time correlation function on the 1-torus. The usual 
time correlation function is defined via 

(xx*) = f dxp(x)xf*(x), (8) 

Jo 

where the time evolution of x is given by the map f(x), 
but here we introduce the time correlation function G(t ) 

G{t)= f dx°p(x 0 )e 2 ^ x °- xt \ (9) 

Jo 

where 

x* = (a*x°) mod 1, (10) 

which respects the fact that the variable x is an angular 
variable on a torus ||. 

For p(x) = 1 we obtain from eqs. (I© 


If the parameters a = m + 2n and e = 2n/(m + 2n) are 
such that both (1 — e)a and ae/2 take integer values m 
and n, then the equation of motion for the coupled map 
system can be written in the compact form 

y; A t jX*\ mod 1, (13) 

where the coupling matrix A has integer elements 

Aij = mSij + n(5i, i+ i + 5 iti - 1 ). (14) 

We will now free ourselves from the specific form jb| ) 
for Ajj , which was physically motivated by the nearest 
neighbor lattice model and show that eq. (0) can 
be solved for any matrix A , which has integer valued el¬ 
ements A^. In order to see this we write © in vector 
notation as 

x t+1 = (Ax*) mod 1, (15) 

where x * = (x\, - ■ ■ ,x* N ) and the modulo is taken for 
each component of the vector Ax*. Then we obtain by 
iterating from the initial condition: 

x t+1 = (dlx*) mod 1 = Ax* — k*, (16) 

where k* is a vector with integer components which rep¬ 
resents the action of the modulo function. This yields 

x‘+ 2 = (Ax t+1 ) mod 1 = (A [Ax* - k‘] ) mod 1 

= (AAx* — A\s*) mod 1 = (A 2 x*) mod 1, (17) 

where the last equality sign only holds because all ele¬ 
ments of the matrix A are integers, such that A\t* is a 
vector with integer components which can be dropped 
under the last modulo function. Since ( 0 ) holds for any 
t, we obtain the closed-form solution as a function of the 
initial value 


c ^ 

G(t)= / dx P (x)e 2 < 1 ~ at ) x = 5 ha t. (11) 

Jo 

In the following section we will demonstrate what 
changes have to be made in order to compute in a similar 
fashion as above the solution of the dynamical equations, 
the invariant density and the time correlation function for 
our coupled map lattice. 


III. A LATTICE OF COUPLED BERNOULLI 
SHIFT MAPS. 


For the Bernoulli shift map /(x) = (ax) mod 1 the 
time evolution for the variables x£ of the coupled map 
lattice becomes according to eq. (G) 


j+i 



( 12 ) 


x* = (A‘x°) mod 1. (18) 

Eq. (0) shows that we can solve not only our cou¬ 
pled map lattice problem ©, but all linearly coupled 
systems, where the coupling occurs via a matrix A with 
integer elements Ajj and the nonlinearity is provided by 
the modulo function. The solution can be obtained by 
first solving the linear problem, i.e. by obtaining A*x° 
and then taking the modulo, which is the same as having 
the linear map acting on an TV-torus in analogy to the 
famous Arnold’s cat map in two dimensions Q . 

Next we investigate the invariant density and the spa- 
tiotemporal correlation function of the coupled map lat¬ 
tice. The first quantity gives us information about the 
measurable time averaged spatial structures in the sys¬ 
tem and the second one tells us about the measurable 
spatiotemporal structures. 
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A. The Invariant Density 


The invariant density p(x) yields the distribution of 
points on the attractor generated by the map x t+1 = 
(Ax*) mod 1. By starting from an initial distribution 
p°(x) it could be obtained as the infinite-time limit of 
/3*(x) in the Frobenius-Perron equation 

p*(x) = J dx'S [x — (A*x') m od l] p 0 (x') (19) 

Since all quantities involved in eq. @ are periodic on 
an IV-torus, the Fourier decomposition of p*(x) contains 
only wavevectors k with integer components, i.e. 

p*(x)=^p*(k)e 2 - kx (20) 

k 

By using the equality exp(27rik ■ [(A*x)mod 1]) = 
exp(27rz[(A t ) T k] • x) eq. ( |2(]| ) yields 

p*(k)=p°((A t ) T k). (21) 

If the initial distribution p°(x) is non-singular, all Fourier 
coefficients vanish for large values of the wavevector, i.e. 

lim p°(k) = 0. (22) 

|k|—>-oo 

For a completely expanding map, where all eigenvalues 
of the matrix A have an absolute value larger than one, 
limt -,00 (A‘) T k = oo for each k ^ 0 and the only non¬ 
vanishing Fourier component becomes 

A°) = J dxp°(x.) = 1, (23) 

which yields a constant invariant density 

p(x) = 1. (24) 

This result is completely analogous to the single map 
case. However we may obtain different results for the 
invariant density if there are contracting directions in the 
phase space. 

Indexing the stable and unstable eigenvalues A and 
right (left) eigenvectors e (e) of the matrix A T with in¬ 
dices “s” and “u” respectively, we have 

(A‘) T k = £ A*(e s ■ k)e s + ]T A* (e“ • k)e“. (25) 

S U 

According to the above, we will only obtain results which 
differ from the trivial expanding case, if there exists at 
least one k / 0, such that its components along the un¬ 
stable directions are all zero, i.e. it is contained in the 
stable manifold W s of the fixed point k = 0 of the “con¬ 
jugate” map 



ki 


FIG. 1. The cat map: none of the integer-component 
wavevectors lies on the stable manifold W s of the fixed point 
k = 0. 

In fact the above argument does not take into account 
one particular feature of the system. Specifically, the 
map ( |26| ) might not be hyperbolic, i.e. it might possess 
the central manifold, defined by the eigenvectors corre¬ 
sponding to |A| = 1. If this is the case (and it is for 
the coupled map lattice © as we shall see below), we 
assume that the invariant density, that we calculate cor¬ 
responds to the physical (or Kolmogorov) measure |j]. 
The latter is calculated by introducing a small amount 
of noise into the system and then taking the zero-noise 
limit. 

Assuming that the physical measure is unique, we can 
obtain the invariant density starting from p°(x) = 1, 
which is equivalent to /3°(k) = 4k.o- Then only the sta¬ 
ble manifold contributes to the invariant density and the 
central manifold can be treated as the unstable one. This 
essentially means, that the averaging over the central 
directions is accomplished by the infinitely small noise 
present in the system. 

On the other hand, eqs. (|l|,^5|) tell us that it is not 
enough to have contracting eigenvalues in order to get a 
non-constant invariant density. 

Let us first consider the case with a single stable di¬ 
rection e s . Since all components ki of a vector k are in¬ 
teger, it is contained within the stable manifold W s only 
if ^k = e s . This in turn means that the components 
{ef, ■ ■ •, e s N } should be mutually rational, i.e. 

e\ : e s 2 : ■ ■ ■ : e s N = k\ : k 2 : • • • : fcjv- (27) 

An example where we have one contracting and one 
expanding direction is the cat map 

A=(; *). (28) 

Although the eigenvalue corresponding to the contracting 
direction is A s = (3 — v / 5)/2 < 1, this map still has 


k t+1 = A t k*. 


( 26 ) 
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a constant invariant density because the components of 
the eigenvector e s = (2,1 — \/5) belonging to A s have a 
non-rational ratio (see fig. [l]) leading to 

lim p*( k) = 5 k ,o- (29) 

t—> OO 


Generally, in order to get a non-constant invariant den¬ 
sity our model must possess a stable manifold, which in 
turn should contain at least one vector with mutually 
rational components. Every vector k* with integer com¬ 
ponents, which is pulled in the long time limit into the 
origin, according to eq. (p5|), 

lim p‘(k*) = lim p° ((AY k*') = p°(0) = 1 (30) 

can be represented as a linear combination of a (usually 
small) number of basis vectors P,j = 1, • • •, M with in¬ 
teger coefficients rij, i.e. k* = J2j n j f J '- The invariant 
density will only contain non-vanishing Fourier compo¬ 
nents with such wavevectors k*, all with weight 1: 


d( x ) = Z 1 


,2 Tri k* 


k* 

M 


M 

=nz* 

3 =1 n 3 


2-Kinn f 3 • 


3 =1 rij 


v 3 / ’ 


(31) 


B. The Coupled Map Lattice 


fZ 2 = (0,-1,0,1, - - -,1); 

f^ /2 = (1,0,-1,0, • • -,0); 

Y /3 = (- 1 ,- 1 , 2 , •■•, 2 ); 

f 2 27r/3 = (-l, 2,-1, - --1); 

P = ( 1,1, - - -, 1). (34) 

Rationality of cos (q*) in not an unexpected result, 
e.g. choosing cos (q*) = —m/2n results in the eigenvalue 
\ q » =0, according to (j32[), which requires 

(f • x 4 ) mod 1 = 0, Vi > 0, (35) 

where we defined f = n q *e q with n q * = N if q* = 
0, 7t/ 2, 7r and 2 N otherwise. This in turn, requires p(x) ~ 
i 5((f • x) mod 1), which is seen to be the case by rewriting 
(|3l|) as 


M 

p (*)=n <5((f J • x) mod 1). (36) 

3 =1 

It is useful to define the projection of the invariant 
density p{x) on a chosen direction g: 

Ps(s) = J S(s- g ■ x)p(x) dx. (37) 

For example, if gi = Sij, eq. (^) gives the distribution 
of the j -th map variable p(xj ) = 1. 

If g coincides with one of the basis directions, i.e. g = 
for some l , the projection 


Up to now our conclusions have been completely gen¬ 
eral for any coupling matrix A with integer elements. 
Let us now consider the condition (^) in more detail for 
our 1-dimensional nearest neighbor model (12). The cor¬ 
responding matrix (0) can be diagonalized by Fourier 
transformation in the space variables i, leading for peri¬ 
odic boundary conditions to eigenvalues 


M 


Pe(s) = IS( s ~vf l -x) F ' x ) dx 

lN j=i p 

= '^2 D p5( s ~i'p), (38) 


(. I N denotes the unit ./V-dimensional cube) becomes sin¬ 
gular: we get a series of 5-functions with an envelope 


A g = to + 2ncos(<7) (32) 

and the corresponding eigenvectors 

e 9 = N^ 1 (cos(q), cos(2 q), • • • , cos (Nq)), 

e q s = 7V _1 (sin(g),sin(2g), • • ■ ,sin (Nq)), (33) 

where q = 2np/N and p = 0, • • •, N/2. 

Of these only a few have mutually rational compo¬ 
nents. For instance, both cos(g) : 1 and sin(2g) : sin(g) 
are rational only if cos (q) is rational, which immediately 
restricts the allowed wavevectors q = 2np/N to a set of 
5 values: q* = 0,7r/3, 7 t/2, 27t/3, n. Each q* generates 
several basis vectors, provided |A g * | < 1: 

f° = (i, 

fJ r/3 = (l,-l,-2,-l,l,2,.--,2); 
f* 73 = (-1,-2,-1,1,2,1, •••,!); 


D p = S(p — i 1 ■ x) TT S(n — • x) dx. (39) 

JlN fin 

Otherwise, the projection (^?[) is a continuous, non¬ 
singular function of parameter s. In other words, only 
the projection on the directions defined by the basis vec¬ 
tors f J is singular. 

In particular, the eigenvector e 9 defines a basis direc¬ 
tion f J if and only if the projection (^) on this eigenvec¬ 
tor (we define p q {s) = p g (s) for g = e 9 ), 

r M 

p q (s)= S(s — e 9 • x) J^[ ^ S(pj — f- 7 • x) dx, (40) 

JlN ,=i Z 

is singular. This implies that e 9 = v q P for some j. 

One can trivially verify that the projection (|4(]) has 
the average 
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FIG. 2. Projection of the invariant density p q (s — s q ): (a) for q ^ q* , arbitrary \ q and also for q = q* , |A ? >1 and (b) for 
q = q*, |A 9 | < 1. We used N = 32. 


s q = J s p q (s)ds = - t S qt o (41) 

and the dispersion given by 

Vq = J ( S ~ s qf P q ( S ) ds = (1 + $q, 0 + <5q,ir) ( 42 ) 


for all q ^ q* and almost always for q = q*. A few 
degenerate cases like p(x) = 8(x i — X 2 ) or p(x) = S(x 1 + 
X 2 — 1) for N = 2, or p(x) = 8(x 1 — 23 ) i5 (£ 2 — £ 4 ) for 
A” = 4 give different dispersions. 

As expected, numerically calculating the projec¬ 
tion p q (s ) on the stable and unstable directions (33), 
we only get a singular distribution for q = q* = 
0, 7t/3, 7 t/ 2 , 27t/3, 7t (fig. |](b)), provided that the respec¬ 
tive eigenvector is stable (|A g | < 1). Otherwise a smooth 
Gaussian-like distribution is obtained (fig. |](a)). 

Indeed one can easily see that in the large-length limit 
both the continuous distribution and the envelope of the 
singular distribution (ptl) become Gaussian: 


p q (s) 


±ch(^) 

( a 1 J 

ii3 i : 


if Vj, e q x f J 7 ^ 0, 


e q = VgV, 


(43) 


where <j>(t) = (2i r) 1 / 2 exp(— t 2 /2) is the normalized 


Gaussian and v q = n q 1 . 


C. The Spatiotemporal Correlations 

The standard spatial correlation function is trivially 
calculated to yield 

C(r) = (xiX i+r ) - (xi)(x i+r ) = S r: 0 (44) 


for the completely expanding case with p(x) = 1 (here 
(•) denotes the average taken with p(x)). 

If there are contracting directions, we rewrite (|d|) as 


iqr 


c (r) = + a lq) e 


= e«,+<, - TLy ,r +4 ^. 0 . («> 


where a s a = a c a = a 0 for all q except a s n = a s n = 0. 
Since a 2 = (1 + 5 qfi + 5 q ^)/24N for all q 7 ^ q *, 

C'(r) = ^ 5r,o + K 2 - ^v) + - 4v)(-l) r 

+ 2 + a c, q - 7777 ) cos(gr). (46) 


3 ’ 2 ’ 3 


This reduces to a ^-correlation (which coincides with 
the result Q) obtained for p(x) = 1 ) in all but a few spe¬ 
cial cases, when a 2 7 ^ (1 +S q fi +S q ^)/24N. For instance, 
choosing m = 0 and n = ±1 yields for N = 4 


p(x) = S(x 1 - £ 3)^(£2 ^ £ 4 ) 


(47) 


and therefore er S 7r / 2 = o p Cj7r / 2 = 0 and a 2 0 = a 2 ^ = 1/24, 
resulting in 


C{r) = — + = { 12 lfr -°’ 2 , 

{ 1 24 24 \ 0 if r = 1,3. 


(48) 


Since the invariant density, although being nontrivial, 
does not tell us much about the spatiotemporal struc¬ 
tures in the system, next we introduce a spatiotemporal 
correlation function Gi(r,t ), which is a straightforward 
generalization of the time correlation function (|jfj): 


G x {r,t)= J dx°p(x°)e 2ni ( x °- x *+r\ 


(49) 
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By expanding p(x) into Fourier series we obtain in 
analogy to (0): 


N 

Gi(r, t)=£ p(k) d(% - + Sij). (50) 

k .7=1 

Since only the non-vanishing Fourier components 
p(k*) = 1 (where k* = n A l ) of the invariant den¬ 
sity (0 contribute, (|50| ) reduces to 

n / m \ 

G i(U *) = II d ( 51 ) • ( 51 ) 

ni,--,nM j=l \Z=1 / 


In a translationally invariant system G^ (r, t) does not de¬ 
pend on i, so we drop the index and fix i (set * = 1 to be 
specific). 

It can be easily verified that the correlation ([49]) is 
short-ranged in both space and time. First we note that it 
vanishes if the vector k(, with components kj = A\ +r ■ — 
<5ij^k>es not lie on the stable manifold W s . According 

A\ +r ,j = Ai ^e} +r e] + ^ e i+r e j ”1-^ > (52) 

where Ai is the largest and A 2 — the next largest eigen¬ 
value and e 1 and e 2 are the respective eigenvectors. For 
increasing t the vector k 4 asymptotically approaches the 
direction defined by e 1 and therefore cannot lie on the 
stable manifold for t > r, where r is some finite (and 
typically small) integer. 

On the other hand for t = 0we have 


G(r,0) 


N 

E IE 

j =1 


, n lfj ^ 1 +r,j 



(53) 


Since all basis vectors ( |34| ) are periodic in space with 
periods 1, 2, 3, 4 or 6 , any linear combination of these 
will also be periodic with period of at most 12. Since the 
vector with components kj = 6i+ r ,j — Si t j is not periodic 
for r 7 ^ 0 , the maximal size of the system with non¬ 
trivial correlation is limited by N = 12. Again, choosing 
m = 0 and n = ±1 for A = 4 as an example, we have 
f 1 = ( 0 ,— 1 , 0 , 1 ) and f 2 = ( 1 , 0 ,— 1 , 0 ) as basis vectors 
and consequently 


G(r,0) 


1 if r = 0,2, 
0 if r = 1 ,3, 


i.e. we retrieve the result ([48]). 


(54) 


IV. DISCUSSION 

To summarize, we have shown that the solution for the 
dynamical behavior of a lattice of Bernoulli maps that 


are coupled by a matrix A with integer coefficients can 
be given in the closed form as x 4 = (A 4 x°)mod 1: the 
dynamical behavior of the coupled map system can be 
described by the repeated action of a linear map A 4 x° on 
variables that are confined to an TV-torus. This picture 
explains that the essentials of the dynamical behavior are 
dictated by the eigenvalues and eigenvectors of A. 

The invariant density p(x) of this system displays 
Fourier coefficients that are different from zero, i.e. is 
non-constant, whenever the stable manifold of the zero 
wavevector contains a non-empty basis of directions f J 
with mutually rational components, generating the infi¬ 
nite asymptotically contracting set of wavevectors. For 
nearest neighbor couplings in a 1-dimensional lattice 
(given by eq. |T~i| ) the maximal number of basis vectors is 
eight (actually even less, since, e.g. |Ao — A^| = 4|n| > 2). 

We have calculated the standard spatial correlation 
function C(r) for the model with nearest neighbor cou¬ 
plings and shown that it is given by C(r) = <5 ri o/12 al¬ 
most always. A few special cases exist however for suf¬ 
ficiently small lattices, where the spatial correlations are 
different. Nevertheless, C(r) always vanishes at suffi¬ 
ciently large distances. 

The invariant density of this system and the spatial 
correlation function display little structure as compared 
to the Lyapunov spectrum, which is, for the nearest 
neighbor coupling, given by A q = log|m + 2ncos(q)\. 
This result shows that the time averaged spatial behavior 
is not simply a straightforward reflection of the Lyapunov 
spectrum (see related work listed in Jl(J). 

We have also calculated the measurable spatiotemporal 
correlation function G(r, t ) for the translationally invari¬ 
ant model and shown that it too is short-ranged in both 
space and time. 

It is instructive to compare our results with the general 
results obtained by Bunimovich and Sinai ||, who proved 
that, for sufficiently small coupling (in our case deter¬ 
mined by parameter e), certain expanding coupled map 
systems with finite-range coupling possess an absolutely 
continuous invariant measure /r(x) : d/x(x) = p(x) dx, 
and also that the time and space correlation functions 
decay exponentially (not slower than exponentially, to 
be exact). 

Our results indicate, that for larger coupling, the in¬ 
variant measure still exists, but might not be absolutely 
continuous due to the fact, that large coupling often 
causes the appearance of contracting directions, even if 
the isolated local maps f(x) are expanding. The space 
and time correlations in our model are seen to decay even 
faster then exponentially, but the few special cases, giving 
non-trivial correlations imply that there might be some 
general relationship between the continuity of the invari¬ 
ant measure and the appearance of coherent structures 
in the system. 

Let us finally point out several directions of further 
research. 

One open problem is the extension of our results to 
higher dimensions and to couplings which have a longer 
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range. In the 1-dimensional case the eigenvectors remain 
also valid for longer-ranged couplings, only the eigenval¬ 
ues change. This means that a model with long, but 
finite, range will have no more structure in the invariant 
density than the short-ranged model. This is of course a 
peculiarity of the Bernoulli shift map, but should again 
be taken as a warning for making conclusions from the 
spatial range of the coupling onto the observable spatial 
patterns. 

Although our solution for the dynamics and the corre¬ 
lation functions hold for general dimensions it would be 
interesting to see what the restrictions on the wavevec- 
tors that generate the basis of the invariant density look 
like in two and three dimensions. 

Finally, one could investigate the dynamical behavior 
of a system, whose time dependence is given a priori by 
the equation ( [Tsl ) also for matrices A with non-integer 
elements. By doing so one will loose the property of the 
original map, that the relation relation (|15j) holds step 
by step, but the trajectories generated by equation (ED 
are well defined. 
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